clear
load('1500.mat')

deltaV = V(end) - Vf*Vs;
deltaH = alt(end) - hf/1000;
deltaS = acos(sin(lat(end)*pi/180)*sin(latf) + cos(lat(end)*pi/180)*cos(latf)*cos(lonf-lon(end)*pi/180))*R0/1000;
deltat = t(end) - tf*ts;